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Strongly coupled QED is a model whose physics is dominated by short-ranged effects. In order to assess which 
features of numerical simulations of the chiral phase transition are universal and which are not, we have formulated 
a quenched version of the model in which photon degrees of freedom are defined on a lattice of spacing a, but 
fermions only on a lattice of spacing 2a. The fermi-photon interaction is then obtained via a blocking procedure, 
whose parameters allow a degree of control over the relative importance of short wavelength modes. Results 
from a variety of models are presented; the critical exponents S and j3 m governing the transition appear to be 
independent of the blocking, or even of whether a gauge- invariant action is used for the photons. 
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In this talk I will describe work done in col- 
laboration with John Kogut. We have been ex- 
perimenting with different variants of the non- 
compact lattice QED action. The motivation is 
twofold: firstly we wish to get a better under- 
standing of chiral symmetry breaking, which oc- 
curs at strong coupling, and in particular to dis- 
entangle the role of various features, eg. whether 
the physics is dominated by short or long ranged 
effects; the relative importance of photon ex- 
change and four-fermi contact interactions; and 
the role, if any, of magnetic monopole excita- 
tions. Secondly, if there were a UV fixed point in 
strongly coupled QED, then there should be the 
possibility of some renormalization group "im- 
provement" which would widen the scaling region 
seen in simulations and give cleaner signals. 

A specific example is the hypothesis that the 
fixed point can be described by the gauged 
Nanrbu Jona-Lasinio model: 
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Studies using the quenched ladder approximation 
|f] |2j predict a line of fixed points running from 
the NJL mean field transition at (a, G) — (0,4) 
to the essential singularity discussed by Miransky 
and Bardeen et al at (a c , 1) with a c = 7r/3. Criti- 
cal exponents of the chiral transition are analytic 
functions of a/a c along the fixed line. Numer- 
ical simulations of quenched non-compact QED 



(NC QED) yield estimates of exponents consis- 
tent with crossing the fixed line at (0.44cv c , 3.06) 
M . Why might the gauged NJL model be a good 
effective theory for lattice QED? One clue is the 
lattice photon propagator. In Feynman gauge, 



D~l{k) 



^5>in 2 ^. 



(2) 



This falls below the continuum form k 2 at the 
edge of the Brillouin zone, and generates an effec- 
tive contact interaction between vector currents, 
as well as the term due to photon exchange. 

Our proposal is to modify the relative weights 
of large and small k modes in the propagator 
by modifying the lattice photon action. The ap- 
proach is simple: from photon fields O^ix) gener- 
ated on a (2A*) 4 lattice using the standard non- 
compact action we form blocked links © M (x) on a 
A^ 4 lattice via: 
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The action remains Gaussian under blocking 
and it is straightforward to examine the new dis- 
persion curves (fig. I). From the large family of 
possible blocked actions we chose two forms to 
examine in simulations: (I) a = 1, b = c = d = 0. 
We found this to enhance short wavelength modes 
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Figure 1. Inverse photon propagator parallel to 
(0,0,0,1) 



over long, effectively increasing the "screening" 
in the model; and (II) a = 1, b = |, c = d = j^. 
This choice gave a relative suppression of short 
wavelength modes. For comparison we also tried 
a "totally screened" random link (RL) model in 
which the {0} are generated independently on 
each link with Gaussian weight exp(— 166^). 

To examine the critical properties of the chi- 
ral transition, we estimated the chiral conden- 
sate (iptp) m the chiral limit bare mass m — > 
on (blocked) lattice sizes 8 4 , 12 4 , and 16 4 in two 
complementary ways pi. First we measure (ipip) 
as a function of m by inverting (^Z>[0] +rn), where 
jp is the staggered fermion covariant derivative, 
using a conjugate gradient alogorithm. We used a 
sequence of m values 0.005, 0.004, . . . , 0.001. Sec- 
ondly, we estimate the spectral density function 
p(X) by using the Lanczos algorithm to find the 
eigenvalues of p[0]- The spectral density is re- 
lated to a branch cut in (?pip(m)} along the imag- 
inary m axis: 

2-k P {\) = Hm(^{-iX + e))-{ijij(-iX-e)). (3) 

On the assumption of power-law scaling at criti- 
cality, we have a relation between the critical cou- 
pling (3 C = 1/ej? and the critical exponent 8: 

(#("»)}| Me =4m 1/S . (4) 

The two parameters can be estimated by looking 



for the best fit of eqn. (4) as a function of (3 to 
the conjugate gradient data. Alternatively, the 
Lanczos data can be fitted using the ansatz 
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As a cross check, using (3) we have the relations 
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By demanding consistency between the two sets 
of measurements, we can achieve accuracy of bet- 
ter than a percent in estimating /3 C , and 5% in 
5. It is also possible to estimate the magnetic ex- 
ponent [3 m from the p(A) data, though with less 
precision: 



p(X = 0) = C((3 C - fff« 
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Another observable we kept track of was the 
size of the largest connected monopole cluster, 
along with the associated cluster susceptibility 
H . By monitoring this on a series of lattice sizes 
8 4 , 10 4 , . . . , 20 4 we were able to estimate the criti- 
cal coupling for the percolation transition. This is 
worth studying, because the percolation thresh- 
old is close to the chiral transition in NC QED, 
and actually appears to coincide once dynamical 
fcrmions are introduced. Although the relevance 
of monopoles to the chiral transition is contro- 
versial [H, the percolation threshold is at least 
a characteristic property of the "QED vacuum" , 
and it will prove interesting in this study. 

Our results for the chiral transition are pre- 
sented for each of the four quenched models we 
have studied in the following table. The suscepti- 
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bility exponent 7 was calculated using the scaling 
relation 7 = /? m (<5 — 1): it is included because the 
quenched ladder approximation predicts it to be 
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Figure 2. Critical inverse propagator /3 C D *(k) 



unity |1| |2| . We see that the main trend is for the 
models which are screened relative to the stan- 
dard formulation (ie. short wavelength modes less 
suppressed) to break chiral symmetry at weaker 
bare coupling, whereas the less screened model 
(Block II) has the transition at stronger coupling. 
The critical inverse photon propagators D~^(k) 
are shown for each model in fig. 2, averaged over 
the corners of the Brillouin zone, since the block- 
ing violates hypercubic symmetry to an extent. 
Also shown for comparison is the critical inverse 
propagator for 4 flavor QED, assuming a critical 
coupling f3 c ~ 0.205 and logarithmic screening 
using the perturbative formula. The crucial ob- 
servation is that all the curves differ strongly in 
the k — > limit, and approach each other near 
the zone edge. This is conclusive evidence that 
chiral symmetry breaking in QED is due to short- 
ranged interactions. It is also striking how sim- 
ilar the estimates for the exponent S are in the 
four models, even for the RL model which vio- 
lates gauge invariance. The universality of 5 and 
(less convincingly) (3 m remains to be understood 
- perhaps it is a feature of the staggered fermion 
formulation. It will be extremely valuable to get 
quantitative results using other formulations. 

In the next table we compare the critical cou- 
plings for the chiral and monopole percolation 
thresholds. From this we learn that the two tran- 
sitions can be made distinct if screening is less- 
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ened: however for "screened" models (Block I, 
Nf ~ 2, 4 QED) the two transitions appear to co- 
incide H II . Figure 2 also supports a crude divi- 
sion of models into "screened" and "unscreened" . 
In the quenched theory, we also see that 



/3 c (xSB) > [3 C (percolation) 
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for all models in which the photon fields are gen- 
erated using a gauge-invariant action. 

The results of this study have raised more ques- 
tions than have been answered. Little support has 
been found for the scenario predicted by the lad- 
der approach to the gauged NJL model: we have 
been unable to affect the critical exponents by 
tuning the short-ranged properties of the model. 
Perhaps a more radical approach, such as a for- 
mulation in momentum space H, is needed. We 
have also shown that monopole percolation can 
be made separate from the chiral transition, at 
least in the quenched case, but have found little 
evidence for the percolation threshold occurring 
at a weaker coupling than the chiral transition, 
except in the RL model. It would appear de- 
sirable to continue to explore this problem from 
many different directions. 
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